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(57) Abstract: A pulse oximetry method and system for improved motion correction provides for the use of a detector output signal 
to obtain a different plurality of differential absorption data sets in corresponding relation to each of a succession of measurement, 
wherein each of the data sets includes differential absorption values for light of a first wavelength and light of a second wavelength. 
The data sets are processed to obtain a relative motion estimate value (126) for each measurement. When the relative motion estimate 
value for a given measurement falls within a predetermined range, a corresponding blood analyte indicator value is adjusted (128) 
in a predetermined manner, wherein the corresponding adjusted blood analyte indicator is employable to obtain at least one blood 
analyte concentration value (130). 
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PULSE OXIMETRY METHOD AND SYSTEM 
WITH IMPROVED MOTION CORRECTION 

FIELD OF THE INVENTION 

5 

The present invention generally pertains to patient monitoring using 
photoplethysmographic devices to generate blood analyte information. More 
particularly, the invention is directed to an efficient and effective approach for 
reducing the undesired effects of motion-contaminated data in pulse oximetry devices. 

10 

BACKGROUND OF THE INVENTION 

In the field of photoplethysmography light signals corresponding with two or 
more different centered wavelengths may be employed to non-invasively determine 

1 5 various blood analyte concentrations. By way of primary example, blood oxygen 

saturation (Sp02) levels of a patient's arterial blood are monitored in pulse oximeters 
by measuring the absorption of oxyhemoglobin and reduced hemoglobin using red 
and infrared light signals. The measured absorption data allows for the calculation of 
the relative concentrations of reduced hemoglobin and oxyhemoglobin, and therefore 

20 Sp02 levels, since reduced hemoglobin absorbs more light than oxyhemoglobin in the 
red band and oxyhemoglobin absorbs more light than reduced hemoglobin in the 
infrared band, and since the absorption relationship of the two analytes in the red and 
infrared bands is known. 

To obtain absorption data, pulse oximeters comprise a probe that is releaseably 

25 attached to a patient's appendage (e.g., finger, ear lobe or the nasal septum). The 
probe directs red and infrared light signals through the appendage, or 
tissue-under-test. The light signals are provided by one or more sources which are 
typically disposed in the probe. A portion of the light signals is absorbed by the 
tissue-under-test and the intensity of the light transmitted through the tissue-under-test 

30 is detected, usually by at least one detector that may be also located in the probe. The 
intensity of an output signal from the detector(s) is utilized to compute Sp02 levels, 
most typically via a processor located in a patient monitor interconnected to the probe. 

As will be appreciated, pulse oximeters rely on the time-varying absorption of 
light by a tissue-under-test as it is supplied with pulsating arterial blood. The 

35 tissue-under-test may contain a number of non-pulsatile light absorbers, including 
capillary and venous blood, as well as muscle, connective tissue and bone. 
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Consequently, detector output signals typically contain a large non-pulsatile, or DC, 
component, and a relatively small pulsatile, or AC, component. It is the small 
pulsatile, AC component that provides the time- varying absorption information 
utilized to compute arterial Sp(>2 levels. 
5 In this regard, the red and infrared signal portions of pulse oximeter detector 

output signals each comprise corresponding large DC and relatively small AC 
components. The red and infrared signal portions have an exponential relationship to 
their respective incident intensities at the detector(s). As such, the argument of the 
red and infrared signal portions have a linear relationship and such portions can be 

10 filtered and processed to obtain a ratio of processed red and infrared signal 

components (e.g., comprising their corresponding AC and DC components), from 
which the concentration of oxyhemoglobin and reduced hemoglobin in the arterial 
blood may be determined. See, e.g., U.S. Patent No. 5,934,277. By utilizing 
additional light signals at different corresponding centered wavelengths it is also 

1 5 known that carboxyhemoglobin and methemoglobin concentrations can be 
determined. See, e.g., U.S. Patent No. 5,842,979. 

As noted, the pulsatile, AC component of a pulse oximeter detector output 
signal is relatively small compared to the non-pulsatile DC component. 
Consequently, the accuracy of analyte measurements can be severely impacted by 

20 small amounts of noise. Of particular concern here is noise that contaminates 

absorption data as a result of undesired variations in the path length of light signals as 
they pass through the tissue-under-test. Such variations are most typically caused by 
patient movement of the appendage to which a pulse oximetry probe is attached. 

A number of different approaches have been utilized to reduce the deleterious 

25 effects of patient motion in pulse oximeters. For example, pulse oximeter probes 

have been developed to enhance the physical interface between the probe and tissue- 
under-test, including the development of various clamp type probe configurations and 
secure wrap-type probe configurations. Further, numerous approaches have been 
developed for addressing motion contaminated data through data processing 

30 techniques. While such processing techniques have achieved a degree of success, 
they often entail extensive signal processing requirements, thereby contributing to 
increased device complexity and componentry costs. 
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SUMMARY OF THE INVENTION 

In view of the foregoing, a general objective of the present invention is to 
provide an improved method and system for addressing motion correction in pulse 
5 oximeters. 

More particularly, primary objectives of the present invention are to provide 
for motion correction in pulse oximeters in a manner that is effective and that reduces 
complexity and component requirements relative to data processing-intensive prior art 
devices. 

10 The above objectives and additional advantages are realized in the present 

invention. In this regard, the present inventor has recognized that patient motion can 
be generally classified into a limited number of motion ranges, or bands, for motion 
correction purposes. A first type of motion, which may be referred to as "low 
motion", corresponds with a range of patient motions that do not have a significant 

1 5 effect on the absorption data comprising a detector output signal. A second type of 
motion, which may be referred to as "clinical motion", corresponds with patient 
motion that primarily affects the AC component of a detector output signal in a 
relatively predictable manner across a range of motion severity. For such clinical 
motion, the present inventor has further recognized that a predetermined adjustment 

20 factor may be effectively utilized to correct blood analyte measurements, thereby 
avoiding complex data processing requirements. 

The above-noted recognitions provide the basis for a number of improvements 
to pulse oximetry systems in which a detector provides an output signal indicative of 
light absorption of a tissue-under-test at each of a plurality of different centered light 

25 wavelengths (e.g. centered at red and infrared wavelengths) and which utilize the 
output signal to obtain blood analyte indicator values for each of a succession of 
measurements (e.g., periodic measurements corresponding with partially overlapping 
or non-overlapping measurement data ) during a patient monitoring procedure. The 
inventive method/system provides for the utilization of the detector output signal to 

30 obtain a corresponding relative motion estimate value for each measurement. For 
such purposes, the detector output signal may be processed to yield a different 
plurality of differential absorption data sets for each measurement, wherein each data 
set includes a first differential absorption value for light at the first wavelength dA;^ 
(e.g., infrared light) and a corresponding-in-time second differential absorption value 
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for light at the second wavelength dA x 2 (e.g., infrared light). As will be appreciated, 
the plurality of data sets corresponding with each given measurement are obtained 
over an associated time period, wherein each successive measurement employs a 
different successive plurality of data sets, and wherein the data sets employed for 
5 successive measurements may or may not partially overlap. By way of primary 
example, the differential absorption data sets may be obtained via derivative or 
logarithmic processing of a series of data samples that correspond with a detected 
infrared portion of a detector output signal and a corresponding-in-time series of data 
samples that correspond with a detected red light portion of the detector output signal. 

10 For each of the measurements the method/system provides for a determination 

of whether or not the corresponding relative motion estimate value (RMEV) is within 
a first predetermined range (e.g. corresponding with a predetermined range of clinical 
motion), wherein for measurements having a corresponding RMEV within the first 
predetermined range the corresponding blood analyte indicator value (BAIV) may be 

1 5 adjusted in a predetermined manner. Such adjustment may provide for the use of a 
predetermined adjustment factor (e.g., determined empirically). For measurements 
indicating clinical motion, the adjusted B AIVs may be utilized to obtain a more 
accurate measure of blood analyte concentration (e.g. the relative concentration of 
oxygenated hemoglobin and reduced hemoglobin for Sp0 2 level determination). As 

20 will be appreciated, one or more clinical motion bands, with corresponding 

predetermined RMEV ranges and adjustment factors, may be employed in accordance 
with the present invention. 

As noted, the present inventor has recognized that for a certain range of patient 
motion, referred to as low motion herein, the impact of the motion on blood analyte 

25 measurements is negligible. As such, the inventive method/system may further 

provide for a determination as to whether the relative motion estimate value (RMEV) 
for each of the plurality of measurements is within a second predetermined range 
(e.g., corresponding with low motion). For each measurement having a 
corresponding RMEV within the second predetermined range, the corresponding 

30 blood analyte indicator value may be employed without adjustment to obtain blood 
analyte concentration values. 

The RMEV obtained with respect to each given measurement period may 
entail the computation of a best-fit function for the corresponding plurality of 
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differential absorption data sets. By way of example, the best-fit function may be 
defined as the slope value of a regression line (e.g. determined via least squares or 
linear regression processing) for a "plot" of the differential absorption values for light 
at the first wavelength dAx l versus the corresponding-in-time differential absorption 
5 values for light at the second wavelength dA x T As may he appreciated, the best fit 
function obtained with respect to each measurement period may also be utilized as the 
corresponding BAIV, e.g., the "best-fit RRatio" value as taught in U.S. Patent No. 
5,954,277. Alternatively, a first plurality of differential absorption data sets may be 
computed for use in the computation of the RMEV; and a second plurality of 

10 differential absorption data sets may be separately computed for use in the 

determination of the BAIVs as taught in a co-pending U.S. patent application entitled 
"Method And Apparatus For Determining Pulse Oximetry Differential Values", filed 
contemporaneous herewith and hereby incorporated by reference. 

In conjunction with the obtainment of a best fit function for each given 

1 5 measurement , the inventive method/system may further provide for the performance 
of a statistical analysis of the corresponding plurality of a differential absorption data 
sets in relation to the best fit function to obtain at least one statistical variance value 
indicative of a degree of any associated motion, wherein the statistical variance 
value(s) may be utilized to compute the corresponding RMEV. In primary 

20 embodiments, the statistical analysis may comprise a principal component analysis 

(PC A). That is, for each measurement the plurality of differential absorption data sets 
(i.e., each comprising corresponding-in-time dA y i 1 and dA^ 2 values) maybe 
statistically analyzed in relation to the corresponding best fit function to obtain at least 
one of a first principal component variance value (Vi) and a second principal 

25 component variance value (V2). In this regard, Vi may be computed to represent an 
amplitude-based statistical variance, while V2 may be computed to represent a scatter- 
based statistical variance. The relative motion estimate value for a given 
measurement may be determined utilizing at least one, and preferably both, of the 
corresponding first and second principal component variance values Vi and/or V2. 

30 More particularly, in one embodiment, for each given one of a plurality of 

measurements the inventive method/system may provide for the calculation of a 
corresponding current motion estimate value (CMEV) utilizing the corresponding Vi 
and V2 values, and the identification of a reference motion value (RMV), wherein the 
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corresponding relative motion estimate value (RMEV) is computed utilizing both the 
CMEV and RMV. In one arrangement, the CMEV for each given measurement may 
be calculated as follows: 

CMEV =V y *V 2 . 

5 In turn, the RMEV for each given measurement may be calculated as a ratio of the 
corresponding CMEV and RMV, e.g.: 

RMV 

The RMV may be established on an ongoing, updated basis to be equal to the 
RMEV that corresponds with the lowest amount of motion (e.g., the lowest RMEV) 
1 0 determined with respect to any prior measurement (e.g. preferably corresponding low 
motion) for a given patient monitoring procedure (i.e., a "low motion reference"). At 
the outset of a given measurement procedure and/or where there is otherwise no prior 
low motion reference basis for establishing an RMV, the RMV for the given 
measurement may be established as follows: 

15 RMV = ^- CMEV*K, 

where Vi, V2 and CMEV are as determined with respect to the given current 
measurement and K is a constant. 

As indicated, for measurements having an RMEV within a first predetermined 
range, i.e., corresponding with clinical motion, a predetermined adjustment factor 

20 (PAF) may be utilized to adjust the corresponding blood analyte indicator value 

(BAIV) for use in blood analyte concentration computations. Such adjustment may 
entail the ready application, e.g., by multiplication, of the PAF to the computed 
BAIV. In this regard, the PAF may be empirically set via statistical analysis of 
clinical motion-affected absorption data and corresponding-in-time, non-motion- 

25 affected absorption data obtained via testing of test subject control groups. 

In one embodiment, for measurements having computed B AIVs that exceed a 
predetermined threshold value, the PAF may be scaled in relation to the 
predetermined threshold value. For example, in applications where each BAIV is 
defined by the slope value of a regression line, the PAF for a given BAIV may be 

30 scaled follows: 
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If A<BAIV<B, then: 



Scaled PAF = 1- 



BAIV-A 
B-A 



*(l-PAF) , 



wherein A and B are predetermined constants, and wherein no scaling occurs when : 



In arrangements where the BAIVs are defined by regression line slope values as noted 
10 above, thp PAF may be preferably set between about .5 and .85, and most preferably 
between about .6 and .75 (e.g., about .6875 in one arrangement where A = .9455 and 
B = .65). 

In additional aspects of the present invention, the inventive method/system 
may further comprise additional features to insure appropriate motion correction 

1 5 where there has been a rapid, significant change in the perfusion of a tissue-under- 
test. By way of example, such variations may occur as a result of the application of a 
tourniquet or blood pressure cuff, or as a result of inadvertent contact with a patient 
"pressure point" (e.g., that cuts off arterial blood flow to the tissue under test). To 
address such instances the inventive method/system may provide for the ongoing, 

20 periodic determination of a motion probability factor (MPF), wherein if the MPF 
exceeds a predetermined threshold the reference motion value (RMV) may be 
adjusted using the MPF. 

The MPF may be computed via a statistical analysis of the plurality of 
differential absorption data sets corresponding with each of all given one of a series of 

25 measurements in relation to a corresponding best fit function computed therefore, 
wherein at least one statistical variance value may be obtained for use as an MPF 
value. More particular, the MPF may be computed utilizing of V 2 and/or Vi values 
computed with respect to each current a series of measurements and computed with 
respect to one or more prior measurements. In one approach, an average V 2 value is 

30 determined in relation to each given measurement, wherein each average V2 value is 
calculated by averaging the sum of the V2 value for a given period and the V2 values 
for a predetermined number of immediately precedent measurements. Then, an MPF 



BAIV>B, 
and wherein PAF may be set to 1 when: 



BAIV<A. 
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for each given measurement may be computed utilizing a comparison, or ratio, 
between the average V2 value computed for the current measurement and the lowest 
average V2 value computed with respect to any prior measurement (e.g., current 
average V2 value/lowest average V2 value). 
5 In one arrangement, the system/method may be established so that an RMV 

adjustment may be made when (i) a predetermined number of successive, non-low 
motion measurements have occurred, and (ii) the motion probability factor (MPF) for 
the current measurement period is determined to exceed a predetermined threshold. 
When the predetermined threshold is exceeded, the RMV may be adjusted as follows: 
1 0 adjusted RMV = CMEV + (CMEV-RMV)(MPF)K, 

wherein CMEV and MPF are as determined with respect to the current measurement 
and K is a constant. 

In addition to the implementation of one or more clinical motion ranges and a 
low motion range, embodiments of the present invention may also address patient 

1 5 motion that exceeds clinical motion and may be referred to as "severe motion". In 
such embodiments, the present invention accommodates an approach wherein no 
blood analyte measurement values are output to a user during "severe motion" (e.g., 
when measurements have corresponding RMEVs outside of the predetermined 
clinical motion and low motion RMEV ranges). 

20 Alternatively, to address severe motion embodiments may be employed which 

provide for the use and adjustment of a previously determined blood analyte indicator 
value (BAIV). For example, the BAIV for the most recent non-severe motion 
measurement may be adjusted by an adjustment factor for use in blood analyte 
concentration computations. The adjustment factor may be based upon tracking of 

25 the DC components of the detector output signal portion(s) corresponding with the 
detected light at the first and/or second eentered-wavelength(s) signals. That is, the 
adjustment factor may be computed as the ratio between such DC component(s) 
measurement, and such DC components corresponding with at least a current, given 
severe motion measurement (i.e. a measurement obtained utilizing data obtained 

30 during severe motion. 

In this regard, in one embodiment the detector output signal may be utilized to 
obtain moving average values of total tissue absorption of the AC and DC signal 
components corresponding with each of the first and second light signals (MAV M , 



8 



WO 02/26126 



PCT/US01/41739 



MAV Z2 ), e.g., for red and infrared signal components, wherein such moving average 
values are computed for a predetermined precedent time interval in conjunction with 
each measurement (e.g., an average for a predetermined number of measurements that 
include the current measurement and a number of prior measurements). In the event 
5 severe motion is detected with respect to a given measurement period, e.g. when the 
corresponding RMEV is not within a predetermined low motion range or any 
predetermined clinical motion range, the BAIV corresponding with the most recent, 
non-severe motion measurement may be adjusted by a DC tracking factor (DCTF), as 
follows: 

MAVxz for most recent pre - severe motion measurement 

10 DCTF= MAV * ATr K, 

— lor current severe — motion measurement 

MAV U 

wherein K is a constant; then: 

adjusted BAIV = BAIV * DCTF. 
The adjusted BAIV may then be employed for enhanced blood analyte determinations 
for the given severe motion measurement. As may be appreciated, this inventive 

1 5 aspect may be advantageously utilized in a number of arrangements, including 
embodiments where no clinical motion bands are defined/employed for motion 
correction purposes. For example, this feature may be separately utilized with the 
system/method disclosed in PCT Publication No. WO 98/04903, PCT Application 
No. PCT/CH97/00282, hereby incorporated by reference. 

20 Additional aspects and other combinations of the present invention will be 

apparent to those skilled in the art upon review of the further description that follows. 



BRIEF DESCRIPTION OF THE DRAWINGS 



25 FIG. 1 is a block diagram of one embodiment of a system for implementation 

of the present invention. 

FIG. 2 illustrates an exemplary pulse waveform received on channel X of the 
FIG. 1 embodiment. 

FIG. 3 illustrates an exemplary pulse waveform received on channel Y of the 
30 FIG. 1 embodiment. 
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FIG. 4 is plot of differential absorption data set values dA x , dA y obtained via 
channels X and Y of the Fig. 1 embodiment, and corresponding with low patient 
motion. 

FIG. 5 is a plot illustrating exemplary 45-degree bias lines for data sets 
5 affected by motion. 

FIG. 6 is a plot illustrating first and second principal components in a principal 
component analysis (PCA) of a set of differential absorption values. 

FIG. 7 is a flow diagram of one embodiment of methodology comprising the 
present invention. 

10 FIG. 8 is a flow diagram of one approach for obtaining relative motion 

estimate values (RMEVs) in the embodiment of FIG. 7. 

FIG. 9 is a flow diagram of one approach for correcting blood analyte 
indicator values (BAIVs) in the embodiment of FIG. 7. 



1 5 DETAILED DESCRIPTION 

Referring to FIG. 1, a pulse oximeter system 100 according to one 
embodiment of the present invention is shown. Included in the system 100 are two 

20 light emitters 102 for radiating red light signals and infrared light signals and a 

photodetectors 104. By way of example, the light emitters 102 may be light emitting 
diodes (LEDs) or laser diodes. 

The emitters 102 and photodetector(s) 104 may be incorporated into a probe 
106 adapted for removable attachment to an appendage 100 of a patient, such as a 

25 finger, earlobe, nasal septum, or other tissue, during a monitoring procedure. In use, 
the probe 106 directs the light signals generated by the light emitters 102 onto one 
side of the appendage 100. Photodetector(s) 104 is positioned on the opposite side of 
the appendage 108 to monitor the intensity of light that is transmitted through the 
tissue and produce an output signal responsive thereto. As may be appreciated, the 

30 red and infrared emitters 102 may be driven by a drive 110 within interconnected 
processor unit 120 so that the output signal of detector(s) 104 is multiplexed 
according to a predetermined scheme (e.g., time-division multiplexed, frequency 
division multiplexed, code-division multiplexed, etc.). 

In turn, a demultiplexer 1 12 may be provided to separate the detector output 

35 signal into infrared and red signal portions on two processing channels, e.g. channels 
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X and Y, respectfully. Signal converters 1 14 may be provided in each channel X, Y, 
with current-to-voltage amplifiers, analog-to-digital converters, and other 
componentry for digitizing the red and infrared signal portions for processing by 
processor unit 120. As will be appreciated, demultiplexer 112 may be provided in a 
5 hardware form or implemented in software downstream of converters 1 14 at 
processor unit 120. 

In accordance with the described system embodiment 100, the processor unit 
120 is provided for automatic motion correction in an effective and efficient maimer. 
For such purposes, processor unit 120 may comprise a front-end module 136 that 

10 receives digitized AC and DC detector signal values on both channels X and Y and 
generates a series of composite infrared and composite red signal values, respectively 
(e.g., 30 values each per sec). Such processing may entail the use of preset AC gain 
and DC gain values, preset A/D ground values, and preset light emitter drive levels to 
normalize the AC and DC values. Processor unit 120 may further comprise one or 

15 more data storage buffer(s) 122 to temporarily store data sets comprising the 
composite red and infrared signal values received from module 136, as well as 
additional values computed by computation modules of processor unit 120. In the 
later regard, the data stored in buffer(s) 122 may be accessed for use by a number of 
preprogrammed computation modules comprising the processor unit 120. 

20 In particular, a blood analyte indicator value computation module 124 may 

access the composite infrared and red data signal values in buffer(s) 122 to compute 
differential infrared and red absorption data sets (e.g., 30 sets/sec.) from which blood 
analyte indicator values may be determined for each of a succession of measurements 
(e.g., every Vi sec.) during patient monitoring. Each blood analyte indicator value 

25 (BAIV) may be defined to corresponded with a best fit function determined for a 
plurality of data sets (e.g., 64 sets) utilized for a given measurement. Computed 
BAIVs may be at least temporarily stored at buffer(s) 122 for use in accordance with 
the described embodiment. 

A relative motion estimate value computation module 126 may also access the 

30 composite infrared and red data signal values in buffer(s) 122 to compute differential 
absorption data sets (e.g., 30 sets sec.) from which relative motion estimate values 
may be successively determined for each measurement (e.g., every Vt. sec.) during 
patient monitoring. In this regard, module 126 may utilize a plurality of differential 
absorption data sets (e.g., 64 sets) for a given measurement to determine a best fit 
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function and corresponding motion variance valves Vi, V2 via a principal component 
analysis (PC A). The motion variance values Vi, V2 corresponding with each given 
measurement period provide an indication of the degree of patient motion occurring 
during the period associated with data sets employed for such measurement. 
5 Module 126 may utilize the Vi and/or V2 value(s) to compute a current motion 

estimate value (CMEV) for each measurement, and to establish a reference motion 
value (RMV). The CMEV and RMV may be employed in module 126 to compute a 
reference motion estimate value (RMEV) for each measurement. Computed RMEV 
and RMV values may be at least temporarily stored at buffer(s) 122. 

10 As indicated by FIG. 1 , the blood analyte indicator values computed at module 

124 and relative motion estimate values computed at module 126 may be provided to 
an adjustment module 128 that may adjust, or correct, the blood analyte indicator 
value (BAIV) for each given measurement when the corresponding relative motion 
estimate value (RMEV) is within a first predetermined motion range corresponding 

1 5 with clinical motion. When the RMEV corresponding with a given measurement 
period falls within a second predetermined range corresponding with low patient 
motion (e.g., non-overlapping with the first range), the BAIV for that measurement 
will not be adjusted by module 128. As will be appreciated, the adjusted and 
unadjusted BAIVs corresponding with successive measurements may be provided to 

20 an analyte computation module 130 for computation of blood analyte concentrations. 
The computed concentration levels may be output to user, e.g., via a display 140 
provided at a patient monitor. The monitor may also house the emitter drive 110, 
signal converters 114, processor unit 120, and in alternate embodiments demultiplexer 
112, emitters 102 and/or detector(s) 104. 

25 Referring further to FIG. 1, processor unit 120 may also include a motion 

probability factor computation module 134 that may interface with module 126 and 
buffer(s) 122 to obtain computed RMEV and PC A variance values Vi and/or V 2 . 
Module 134 may employ such values to compute a motion probability factor (MPF) 
on an ongoing updated basis (e.g. corresponding with each measurement). In turn, the 

30 MPF may be utilized to adjust the RMV utilized in module 126 to compute the 
RMEVs. 

Finally, the processor 120 may further include an average absorption 
computation module 132 which may access the composite infrared and red values in 
buffer(s) 122 to compute average total infrared and red absorption values for 
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predetermined intervals(e.g., corresponding with a predetermined number of 
successive measurements), in corresponding relation to each given measurement. As 
will be further described, such average absoiption values may be utilized in certain 
embodiments in connection with the correction of blood analyte indicator values 
5 utilized for measurements having corresponding relative motion estimate values that 
fall within a severe patient motion range. 

As noted, modules 124 and/or 126 may access data buffer(s) 122 to obtain 
composite infrared and red values from which differential absorption data sets may be 
computed with respect to each given measurement. For purposes of describing such 

10 differential absorption data sets, reference is now made to FIGS. 3 and 4 which 
illustrate pulse waveforms 300 and 400, corresponding with the intensity of the 
digitized, corresponding-in-time composite infrared and red signal values generated at 
module 136 of processor unit 120 for channels X and Y, respectively. By way of 
example, derivative processing of successive data samples on channel X and on 

1 5 channel Y may be approximated to obtain differential absorption valves dA x and dA y , 
as follows: 

dA x s DA X = (I t , x - I MjX ) / [(I tjX + I M , X ) /2], and 
dA y = DA y = (I t , y - I t _ a , y ) / [(I t , y + I M , y ) / 2]. 
Alternatively, the differential absorption values may be obtained by calculating 
20 differences between the logarithm values of successive infrared and successive red 
composite signal values. 

In one embodiment, the computation of differential absorption valves 
dA x , dA y at modules 124 and 126 may be conducted separately, wherein module 126 
employs processing techniques taught, in a copending U.S. Patent Application entitled 
25 "Method and Apparatus for Determining Pulse Oximetry Differential Values", filed 
contemporaneous herewith and hereby incorporated by reference. Alternatively, 
modules 124 and 126 may utilize the same dA x , dA y values, with best fit functions 
computed for each measurement utilized both modules 124 and 126. 

To further describe such best fit functions, reference is now made to FIG. 5. 
30 In FIG. 5, a plurality of differential absorption data sets dA x , dA y that correspond with 
low motion are presented. That is, each data point represents a plot of a differential 
absorption value dA x versus a corresponding-in-time differential absorption value 
dA y . A linear regression analysis has been performed on the data set to determine a 
regression line 500, or best fit function, that best fits the data sets. The slope value of 
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line 500 represents a normalized ratio of dA x and dA y for the corresponding 
measurement and may be used in a known manner to determine an Sp02 level. 

As noted, module 126 of processor unit 120 computes a relative motion 
estimate value (RMEV) for each measurement. Each such value may be determined 
5 utilizing a computed best fit function for the corresponding differential absorption 
data sets. For purposes of further explanation, reference is now made to FIG. 6 which 
illustrates in plot form an important observation of how motion may affect the data 
corresponding with a given measurement period. In particular, it has been found that 
for motion within a clinical motion range the effects upon corresponding differential 

10 absorption values dA x and dA y may be substantially the same. Thus a plotted data 
point affected by motion may be biased, e.g., along a 45-degree angle line. The 
regression line 600 shows where plotted data points should lie if no motion is present. 
The lines slanted at 45-degree 604 represent how data points 602 may be displaced 
when the differential absorption values dA x and dA y are affected by clinical motion. 

15 In one embodiment, the operation of module 126 is based upon the 

observation described above in relation to FIG. 6. More particularly, module 126 may 
utilize a best fit function, or regression line, for each given measurement to perform a 
principal component analysis (PCA) of the corresponding differential absorption data 
sets, wherein variance values Vi and V2 may be determined. FIG. 7 graphically 

20 illustrates the results of a PCA performed on a plurality of data sets, having a best fit 
line 700, to obtain corresponding Vi and V 2 values. The distribution box 702 of the 
data points (not shown individually) may be described using first and second 
principal component vectors and derived from the PCA. Each of the principal 
component vectors accounts for a different amount of variation among the data points. 

25 The first principal component vector describes the cluster of data points disposed 
along a longitudinal axis 704. The second principal component vector corresponds 
with an axis 706 extending perpendicular to the first axis 704. With vectors 
describing variations among the set of data points relative to axes 704, 706, any 
number of algorithms may be used to estimate an amount of motion associated 

30 therewith. 

For example, in one approach, let x denote the dA for channel X and y denote 
the dA for channel Y. The calculations may be based on a buffer of n pairs of data 
(xi,yi). The calculations may be based on the five following summations: 
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1=1 

S 2 =£x, 2 , 



s 3 = £y <> 

/=1 



s 4 =£y?, 

1=1 

and 

These are then used to calculate 



u = S 2 - 5l 



10 and 

P = S 5 -^A. 
The regression slope may then calculated to be 



w—u + - w) 2 + 4» i? 2 
1 5 The correlation coefficient between x and y is 



r = P . 

The variation of the points along the regression line is the variance Vi associated with 
20 the first principal component, and the variation of the perpendicular distance to this 

line is the variance V2 associated with the second principal component. As previously 
noted, these values are useful in evaluating the presence and degree of motion. 
The variance Vi of the first principal component may be expressed as: 
v = u + J3 2 w+2/3p 

1 uT 2 ' 
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and the variance V2 of the second principal component may be expressed as: 
_ w+ J3 2 u-2fip 

2 TTp ■ 

The utilization of Vi and V 2 values in conjunction with various modules comprising 
processor 120 will be further described below. 
5 In this regard, FIGS. 7-9 illustrate process steps in the operation of one 

implementation of the present invention. As shown in FIG. 7, a detector output 
current signal may be initially converted (step 202) into a digital form and composite 
red and infrared values may be obtained from the converted signal. Following signal 
conversion, the red and infrared detector output signal portions may be processed to 

10 obtain at least a first plurality of differential absorption data sets for each successive 
measurement (step 204). A plurality of the data sets may be employed to obtain a 
corresponding blood analyte indicator value (BAIV) for each measurement (step 206), 
e.g., a best fit function as described above. Correspondingly, a plurality of data sets 
may be utilized to obtain a relative motion estimate value (RMEV) for each 

1 5 measurement (step 20 8). 

In one embodiment, the RMEV for each measurement may be determined 
according to the method steps shown in FIG. 8. As illustrated therein, a principal 
component analysis (PCA) may be conducted on the data sets corresponding with 
each current measurement relative to a corresponding computed best fit function to 

20 obtain variance values Vi, V 2 (step 220). In turn, a current motion estimate value 
CMEV may be determined with respect to each current measurement (step 222), as 
follows: 

CMEV = V X *V 2 . 

As further shown in FIG. 8, a reference motion value (RMV) may also be established 
25 on an ongoing updated basis. More particularly, the RMV may be established to be 
equal to the lowest RMEV for any prior measurement having an RMEV in the low 
motion range (steps 224, 226). If no prior low motion measurement is available (step 
224), the RMV for a given measurement may be determined as follows (step 228): 

RMV = ^- CMEV K , 

30 where Vi, V2 and CMEV are determined in the manner described above for the 
current measurement, and K is predetermined constant. 
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Pursuant to the determination of Vi and V2 values, the illustrated process may 
also provide for the computation of a motion probability factor (MPF) on an ongoing, 
updated basis utilizing the V2 measurement values (step 230). More particularly, an 
average V2 value may be computed for each measurement by averaging the V2 value 
5 for the given measurement together with the V2 determined with respect to 

predetermined number of immediately precedent measurements. In turn, the MPF at a 
given measurement may be determined by comparing the average V 2 value for the 
current measurement with the lowest average V2 value computed with respect to any 
prior measurement (e.g., by computing a ratio value therebetween). 

1 0 The MPF may be advantageously utilized to adjust the RMV (e.g., as 

determined in step 226) to address rapid changes in the perfusion of the tissue-under- 
test. In the illustrated process embodiment, such adjustment may be made in the 
event that the MPF exceeds a predetermined threshold after a predetermined number 
of successive non-low motion measurements (step 232). In the event an adjustment to 

1 5 an RMV is to be made, the MPF and RMV may be utilized together with the CMEV 
for the current measurement to make such adjustment (step 234). More particularly, 
such adjustment may be made as follows: 

A Adjustment = (CMEV - RMV) * MPF, 

and, 

20 if i A Adjustment < CMEV * A, 



then: 



and, 



A Adjustment = CMEV * A, 



RMV = CMEV - A Adjustment, 
25 wherein A is a predetermined constant. 

As illustrated in FIG. 8, the CMEVs for each of the measurements, together 
with the RMV (as set in Step 226, Step 228, or as adjusted in Step 234) may be 
utilized to compute the RMV for each measurement (step 236). In the illustrated 
embodiment, the RMV values may be determined as follows: 

30 RMEV= ™EL , 

RMV 

wherein CMEV corresponds with the current measurement. 

Returning now to FIG. 7, the illustrated process further provides for a 
determination of whether or not the RMEV for each given measurement is within a 
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predetermined low motion range (step 210). As previously described, the 
predetermined low motion range may be set so that no adjustment of the blood analyte 
indicator value computed at step 206 is necessary (step 212). In the event that the 
RMEV for a given measurement does not fall within the low motion range, the 
5 corresponding blood analyte indicator may be corrected (step 214). For such 

purposes, reference will now be made to FIG. 9, which illustrates one embodiment for 
correction. 

In particular, and as shown in FIG. 9, an initial determination may be made as 
to whether the relative motion estimate value (RMEV) for a given measurement is 

10 within a predetermined clinical motion range (step 240). In the event clinical motion 
is identified, the corresponding blood analyte indicator value (BAIV) may be adjusted 
in a predetermined manner. In one embodiment a best fit function, or slope value, 
computed for the given measurement may be simply adjusted by a predetermined 
factor to obtain an adjusted BAIV (step 242). By way of example, in one 

1 5 embodiment if a given BAIV is greater than a predetermined constant B the BAIV 
may be multiplied by a predetermined adjustment factor (PAF) to obtain the adjusted 
BAIV. On the other hand, if the BAIV is less than the predetermined constant B yet 
greater than a predetermined constant A, the PAF may be scaled as follows: 



20 In the event that the BAIV for a given measurement is less than the predetermined 
value A no adjustment of the BAIV is conducted. 

In the event that a given relative motion estimate value (RMEV) is not within 
the second predetermined range, a number of options are possible. In one 
embodiment, a corrected blood analyte indicator may be obtained utilizing a BAIV 

25 for a prior low motion measurement (step 244), as adjusted by DC tracking 

component (step 244). More particularly, the BAIV computed for the most recent 
non-severe motion measurement preceding a given severe motion period may be 
utilized in the computation of a BAIV. In particular, the pre-severe motion 
measurement BAIV may be multiplied by a DC tracking factor (DCTF) to obtain a 

30 BAIV for the severe motion measurement. The DCTF may be computed as follows: 
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jor most recent pre-severe motion measurement 

DCTF = MAV K ' 

— for current measurement 

MAV n 



Wherein MAV xx , MAV A2 values are average total absorption values for the infrared 
and red portions of the detector output signal, respectively, as computed for the pre- 
5 severe motion measurement and for the current severe motion measurement, and K is 
a predetermined constant. To compute MAV A1 MAV X2 the light absorption values 
(i.e., DC plus AC components) for the infrared and red portions, respectively, of the 
detector output signal may be averaged for a predetermined number of successive 
measurements. 

10 Returning now to FIG. 7, the illustrated embodiment finally provides 

for the use of blood analyte indicators for each of the measurement periods, as either 
adjusted or unadjusted, to compute a blood analyte concentration (step 212 and 216). 
By way of example, Sp0 2 levels may be computed in a manner as taught in U.S. 
Patent No. 5,934,277. In turn, the computed Sp0 2 values may be output to a display 

15 1 40 as previously noted (step 218). 

The embodiment described above is for exemplary purposes only and is not 
intended to limit the scope of the present invention. Various adaptations, 
modifications and extensions of the described system/method will be apparent to 
those skilled in the art and are intended to be within the scope of the invention as 

20 defined by the claims which follow. 
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CLAIMS 

1 . A method for use in a pulse oximetry system which provides a detector 
output indicative of light absorption by a tissue-under-test at each of a plurality of 
different light wavelengths and which employs the detector output to compute blood 
5 analyte indicator values for each of a plurality of measurements, comprising: 

utilizing said detector output to obtain a corresponding relative motion 
estimate value for each of said plurality of measurements; and, 

determining whether the corresponding relative motion estimate value for each 
of said plurality of measurements is within a first predetermined range, wherein for at 
10 least one of said plurality of measurements having a corresponding relative motion 
estimate value within the first predetermined range the corresponding blood analyte 
indicator value is adjusted in a predetermined adjustment manner, and wherein said 
adjusted blood analyte indicator value is employable to obtain a blood analyte 
concentration value. 

15 2. A method as recited in Claim 1 , wherein for said at least one of said 

plurality of measurements having a corresponding relative motion estimate value 
within the first predetermined range said predetermined adjustment manner includes: 

adjusting the corresponding blood analyte indicator value utilizing a 
predetermined adjustment factor. 
20 3. A method as recited in Claim 2, wherein said adjusting step includes: 

multiplying the corresponding blood analyte indicator value by said 
predetermined adjustment factor. 

4. A method as recited in Claim 3, wherein said adjusting step further 
includes: 

25 comparing the corresponding blood analyte indicator value to a predetermined 

threshold value, wherein said predetermined adjustment factor is scaled when said 
blood analyte indicator value exceeds said predetermined threshold value. 

5. A method as recited in Claim 3, wherein said predetermined 
adjustment factor is set to a value between about .5 and .85. 

30 6. A method as recited in Claim 1, said utilizing step including: 

processing said detector output to obtain a different plurality of data sets 
corresponding with each of said plurality of measurements, each of said plurality of 
data sets comprising a first differential absorption value for light at a first wavelength 
and a second differential absorption value for light at a second wavelength, wherein 
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for each of said plurality of measurements a corresponding best fit function is 
computed for the corresponding plurality of data sets. 

7. A method as cited in Claim 6, wherein for each of said plurality of 
measurements said utilizing step further includes: 
5 performing a statistical analysis of the corresponding plurality of data 

sets in relation to the corresponding best fit function to obtain at least one statistical 
variance value indicative of a relative degree of patient motion; and, 

computing said corresponding relative motion estimate value using 
said corresponding at least one statistical variance value. 
10 8. A method as cited in Claim 6, wherein for each of said plurality of 

measurements said utilizing step further includes: 

performing a principal component analysis of the corresponding plurality of 
data sets in relation to the corresponding best fit function to obtain at least one of a 
first principal component variance value (Vi) and a second principal component 
1 5 variance value (V 2 ); and 

computing said corresponding relative motion estimate value using at least one 
of said corresponding first principal component variance value Vi and said 
corresponding second principal component value V2. 

9. A method as recited in Claim 8, wherein for each given one of said 
20 plurality of measurements said utilizing step further includes: 

calculating a corresponding current motion estimate value (CMEV) using at 
least one of said corresponding Vi and V2 values; 

establishing a reference motion value (RMV); and 

computing said corresponding relative motion estimate value (RMEV) using 
25 said corresponding CMEV and said RMV. 

10. A method as recited in Claim 9, wherein for each of said plurality of 
measurements: 

said CMEV is calculated as follows: 

CMEV = Vj*V 2 ; and, 

30 said corresponding RMEV is determined as follows: 

CMEV 
RMEV= rmev ' 
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11. A method as recited in Claim 10 5 wherein for at least one of said 
plurality of measurements establishing step includes: 

setting RMV = CMEV * K 9 

Vi, V 2 and CMEV are as determined for said at least one of said plurality of 
5 measurements, and wherein K is a constant. 

12. A method as recited in Claim 9, wherein for at least one of said 
plurality of measurements said establishing step includes: 

setting the RMV equal to a lowest RMEV corresponding with a prior one of 
said plurality of measurements. 
10 13. A method as recited in Claim 12, wherein said establishing step further 

includes: 

computing a motion probability the using factor (MPF) on an ongoing, 
updated basis utilizing at least one of said Vi values and said V2 values, and, MPF to 
adjust said RMV when the MPF exceeds a predetermined threshold. 
15 14. A method as recited in Claim 13, wherein said computing step further 

includes: 

determining an average V2 value in relation to each given one of a series of 
measurements by averaging a sum of the V2 value corresponding with the given 
measurement and the V2 values corresponding with a predetermined number of prior 

20 measurements comprising said series, and, 

comparing the average V2 value corresponding with each given one of the 
series of measurement periods with the lowest average V2 value corresponding with 
any precedent one of said series of measurement periods to obtain a ratio 
therebetween, wherein said ratio is employable as said MPF. 

25 1 5. A method as recited in Claim 14, wherein in said using step said RMV 

is adjusted to an adjusted RMV as follows: 

adjusted RMV = CMEV + (CMEV-RMV) (MPF) K, wherein CMEV and MPF areas 
determined in relation to a given current measurement and is a predetermined 
constant. 

30 16. A method as recited in Claim 1, further comprising: 

determining whether the relative motion estimate value corresponding with 
each of said plurality of measurements is within a second predetermined range, non- 
overlapping with said first predetermined range, wherein for each measurement 
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having a corresponding relative motion estimate value within said second 
predetermined range the corresponding blood analyte indicator value is employable to 
obtain at least one blood analyte concentration value free from adjustment. 

17. A method for use in pulse oximetry system which provides a detector 

5 output indicative of light absorption by tissue under test at each a plurality of different 
light wavelengths and which employs the detector output to compute blood analyte 
indicator value for each of the plurality of measurements, comprising: 

processing said detector output to obtain a plurality of data sets corresponding 
with each of said plurality of measurements, each of said plurality of data sets 
10 comprising a first differential absorption value for light at a first wavelength and a 
second differential absorption value light at a second wave wavelength, wherein for 
each of said plurality of measurements corresponding best fit function is computed for 
the corresponding plurality of data sets; 

performing for each of said plurality of measurements a statistical analysis of 
15 the corresponding plurality of data sets in relation to the corresponding best fit 
function to obtain at least one statistical variance value; 

computing a relative motion estimate value for each of said plurality of 
measurements using said corresponding at least one statistical variance value; and 

determining whether the corresponding relative motion estimate value for each 
20 of said plurality of said measurements is within a first predetermined range, wherein 
for at least one of a plurality of measurements having a corresponding relative motion 
estimate value within the first predetermined range, the corresponding blood analyte 
indicator value is adjusted utilizing a predetermined adjustment factor, and wherein 
the adjusted blood analyte indicator value is employable to obtain a blood analyte 
25 concentration value. 

18. A method as provided in Claim 1 6, wherein said blood analyte 
indicator values are each defined by a slope value of a regression line computed with 
respect to the corresponding plurality of data sets for a given measurement, and 
wherein the predetermined adjustment factor is set to a value between about .6 and 

30 .75. 

19. A method as provided in Claim 1 7, wherein for each of said plurality 
of measurements said performing step further includes: 

conducting a principal component analysis of the corresponding plurality of 
data sets in relation to the corresponding best fit function to obtain at least one of a 
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first principal component variance value (Vi) and a second principal component 
variance value (V2); and, 

computing said corresponding relative motion estimate value using at least one 
of said corresponding first principal component variance value (Vi) and said 
5 corresponding second principal component variance value (V2). 

20. A method as recited in Claim 19, wherein said method further 
includes: 

establishing a reference motion value (RMV); and, 
wherein for each given one of said plurality of measurements said computing step 
1 0 further includes : 

calculating a corresponding current motion estimate value (CMEV) using at 
least one of said corresponding Vj and V2 values; and, 

computing said corresponding relative motion estimate value (RMEV) using 
said corresponding CMEV and said corresponding RMEV. 
15 21. A method as recited in Claim 20, wherein said establishing step further 

includes: 

computing a motion probability factor (MPF) on an updated basis utilizing at 
least one of said Vi values and said V2 values; and, 

adjusting said RMV utilizing said MPF when said MPF exceeds a 
20 predetermined threshold. 
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